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Abstract 
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^^ , We calculate the two form-factors for the four Kaon to pion transitions via a 

vector current to order p^ in Chiral Perturbation Theory to first order in isospin 
breaking via the quark masses. In addition we derive relations between these form- 
factors valid to first order in the up-down quark-mass difference but to all orders in 
5^ I Chiral Perturbation Theory. 

We present numerical results for all eight form-factors at t = and for varying t 
and for the scalar form-factors at the Callan-Treiman point. 
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1 Introduction 

The seinileptonic decays of a Kaon to a pion and two leptons play a significant role in flavour 
physics. On the one hand, the weak decays Ki^ are a main source of our knowledge of the 
CKM matrix element Vus and on the other hand, the rare decays to a lepton-anti-lepton 
or neutrino-antineutrino pair provide a good testing bed for loop effects in flavour physics. 
The form-factors themselves quantize the hadronic uncertainties as can be exemplified by 
the so-called master formulae. This is e.g. for i^^s, see [Ij and references therein. 
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A similar formula exists for the rare decays, see p] and references therein. We will refer to 
the Kis decays as weak decays and the ones with a lepton-antilepton pair as rare decays. 

Theoretical work on these form-factors goes back a long way. In Chiral Perturbation 
Theory (ChPT), the lowest-order (LO) result dates back to while the next-to-leading- 
order (NLO) was evaluated by Gasser and Leutwyler ^. They calculated the vector 
form-factor /+ for the weak decays including the isospin breaking due to rriu — rriii and the 
scalar form-factor /o in the isospin limit. The form-factors are known in the isospin limit 
to next-to-next-to-leading order (NNLO) in Chiral Perturbation Theory (ChPT) J5i|6j. In 
[5] a comparison with experimental results was done and a useful relation for the order p^ 
constants needed for /+(0) obtained. The isospin breaking to NLO for the vector form- 
factors /+(t) for the rare decays was done in [2| to NLO. The electromagnetic corrections 
to NLO, i.e. order e^p^, are also known, for the weak decays [7| and for the rare decays [2]. 
In this paper we calculate the isospin breaking corrections due to the quark-mass difference 
rUu — TTid to the vector and scalar form-factors to NNLO order in ChPT for all eight form- 
factors. The NNLO results are new for all form-factors while the NLO results are new 
for the scalar form- factors. Some preliminary results were reported in [8j. In addition we 
discuss the results on ratios of form-factors to NNLO. Some of these ratios were observed 
to have special features at NLO in [1] and [2]. We prove that the relations f llip and (1121) 
are valid to all orders in ChPT to first order in m„ — m^. The double ratio flT2|) was also 
discussed in [2] but not proven there. There exists also work using dispersion relation for 
the form-factors in the isospin limit, see [9l [10] and references therein. 

This paper is organized as follows. We define the form-factors and derive the relations 
the form- factors should satisfy to all orders in ChPT and first order in m„ — irid in Sect. |2J 
Next we give a short discussion of ChPT in Sect. [3] and derive how tt^-t] mixing can be 
taken into account to NNLO in ChPT in Sect. HI Sect. |5] defines the various ratios of form- 
factors we use and discusses how they are obeyed at NLO and NNLO. We also discuss 
there some general aspects of our calculation and give the LO results. Explicit formulas 
are not provided at NNLO, they are simply too long but we present the NLO formulas 
in App. |A] and the dependence on the order p^ low-energy constants (LECs) in App. |Bl 
The estimate of these order p^ LECs we use is presented in Sect. [61 Our main results 
are presented numerically in Sect. [71 These include, numerical results on the values of 



/+(0), its t- dependence, ratios as a function of t and the deviation from Fk/ Ft, at the 
Callan-Treiman point. A short summary is given in Sect. [HI 



2 Form-factors and isospin relations 

In this paper we deal with the four matrix-elements 

{7:%p')\-s^,umK+{p)) - 

{n-{p')\-s^,umK\p)) - 
(vr+(p')|^7,^(0)|K+(p)) = 

We have thus in total a set of 8 form-factors. They depend on 

t={p'-p)\ 

the square of the four momentum transfer to the leptons. The form-factors are normalized 
such that 

/f-^(0) = l (7) 

in the SU{?)) limit of m^ = m^ = rUs- In the isospin limit 
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jK-K jg Tj-eieTTed to as the vector form-factor, because it specifies the P-wave projection of 
the crossed channel matrix-elements (s7^g(0)| | K^,tt^ in >. The S-wave projection is 
described by the scalar form-factor 
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We will refer to the decays as the charged weak for ([2]), neutral weak for ([3]), charged rare 
for dH) and neutral rare for ([5]). 

In this paper we derive the isospin breaking due to the quark-mass difference ttIu — Tn^ 
to NNLO for the eight form- factors defined above. We do this to first order in isospin 
breaking. Let us now derive first some general properties. The isospin-breaking operator 
(1/2) (m„ — rriii) (uu — dd) has isospin one. The pions have isospin one and the Kaons as 
well as the vector operator are in an isospin 1/2 multiplet. To first order in isospin breaking 
from 6 = rriu — rrid the form-factors described above can be rewritten in the form 
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for £ = +,—, 0. The form ( ITO|) is a direct consequence of the Wigner-Eckart theorem. This 
can be interpreted as that the size of isospin breaking depends on the final pion and the 
sign also depends on which kaon is in the initial state. 
As a consequence of ( ITOl) we obtain the relations 



fi (t) - fr~{t) - in (t) + fF^ it) = + o{6') , (11) 

and 



rit)^ iKo '%..:::. = ^+o{s^) (12) 



/r-°(t)/r-"(t) 
/r-(t)/r-"(t) 

These relations do not have to be satisfied when electromagnetic corrections are included. 
Photon exchange contains isospin 0, 1 and 2 parts allowing different corrections to all four 
amplitudes. The isospin and 1 parts do satisfy the same relations, but not the isospin 2 
part. 

The relations are valid for all three form-factors ff^''^^ with £ = +,—, 0. They are also 
valid if the currents in ([2115]) are replaced by the scalar densities su{0) and sd{0). 

3 Chiral Perturbation Theory 

Chiral Perturbation Theory (ChPT) is an effective field theory to describe the strong 
interactions at very low energy. The effective Lagrangian is constructed based on two im- 
portant properties of the physical hadron spectrum. Pseudo-scalar mesons, the lowest-lying 
states in the spectrum are separated from the rest of the hadrons, i.e. there exists a mass 
gap. This allows the heavier particles to decouple from the dynamics of the pseudo-scalar 
mesons. Their influence can be described by point-like couplings. The other important fact 
is that the spectrum does not show the chiral symmetry of the underlying theory (QCD). 
The pseudo-scalars are assumed to be the the pseudo-Goldstone particles emerging from 
the spontaneous breaking of this chiral symmetry. The nonzero but small mass of the 
pseudo-scalar mesons are because quarks have a finite mass in, reality which breaks the 
chiral symmetry explicitly. 

According to the Goldstone's theorem, the Goldstone particles do not interact at zero 
momentum. This immediately offers a weakly interacting theory as a basis for perturbation 
theory. The first systematic consideration on the applicability of the effective Lagrangians 
was made by Weinberg [TT] and Gasser and Leutwyler [12] . The effective chiral Lagrangian 
is an expansion in momentum and quark masses. In the chiral power-counting, quark 
masses are of order p'^. Taking into account the Lorentz invariance and chiral symmetry, 
the lowest order chiral Lagrangian which also complies with the discrete symmetries can 
be written down as 

A = ^ {{D^UD^'U^) + {xU^ + Ux^)) (13) 

and the next-to-leading Lagrangian with the introduction of the external field technique 
was written down by Gasser and Leutwyler [I3] and reads 



£4 = Li{D^U^D^U)^ + L2{D^U^DM){D^U^D''U) 



+L,{D>^U^D^U{x'U + U^x)) + L^{X^U + xU^? 
+Lj{x^U - xU^? + L,{x^Ux^U + XU^XU^) 
-iLg{F^^D''UD''U^ + FJ^^D^U^ D^U) 

+L,,{U^F;iUF'^n + H,{F;^^F^^^ + F^^^F^^O + H^{xh) ■ (14) 

The matrix U G SU{3) contains the pseudo-scalars and its exponential representation is 

f/(0) = exp(zv^0/Fo) , (15) 



where 
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The external fields are defined through the covariant derivatives and field strength tensor 
as 

D„U = d„U - ir,U + lUl, , F,i = d^l, - dj^ - t [l^, Q , (17) 
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The right-handed and left-handed external fields are denoted by r^ and l^ respectively. The 
Hermitian 3x3 matrix x contains the scalar (s) and pseudo-scalar (p) external densities and 
is given as x = 2i?o (s + ip). The constants Fq and Bq are related to the pion decay constant 
and quark condensate respectively. There are however, 10-1-2 unknown free parameters in 
the Lagrangian £4 where these effective constants contain the effects of heavy degrees of 
freedom and can be determined by invoking experimental data as well as by Lattice QCD 
technique. One of the theoretical approach, on the other hand, is the application of the 
resonance chiral perturbation which provides an approximate estimate of the low energy 
constants (LECs). The extention of the chiral Lagrangian to the next-to-next-to- leading 
order is also accomplished [H]. At this order there are a large number of LECs, 90-1-4. 

The external scalar field s contains the quark masses and the mass terms in the lowest 
order Lagrangian £2 can be diagonalized exactly. In the presence of rriu 7^ Tna the physical 
tt" and rj differ from the triplet and octet states via a lowest-order mixing angle e as 
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The lowest order mixing angle is 
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A review on ChPT to order p^ is [15] . References to other recent reviews and lectures 
can be found there. 




Figure 1: The full three-point Green function is represented. The oval stands for the 
amputated three-point Green function and circles indicate the full two-point functions. 
The solid lines are external mesons and the dashed line labeled by i , indicates the sum 
over states implied in the external leg where mixing occurs. The wiggly line indicates the 
external vector current. 

4 Matrix-elements in the presence of mixing 

For this work we need to work out the matrix elements defined earlier in the presence 
of mixing. These matrix elements can be determined from three-point Green functions. 
Two of the external legs are the meson propagators and the third one is the external field. 
The matrix element is obtained from the Green function using the Lehmann-Symanzik- 
Zimmermann (LSZ) reduction formula. The matrix element is related to the residue of 
the Green function in momentum space where all the propagators are continued to the 
on-shell mass. The case of two point-function with one leg undergoing mixing is worked 
out in [inj and we generalized this to a four point-function with mixing on two external 
legs in [TTj. In this article we study the form-factors in K^ — > vr-' transitions where, in case 
of the neutral pion in the decay product, mixing should also taken into account. Fig. [1] 
depicts the three-point Green function relevant for this work where we have only considered 
mixing in one external propagator. Amplitudes are obtained via 

/ i-iY' \ " 
A....i. = , ^ ^ n ii"^,(^' - ^') Gi,....„(A:i, . . . , A;„) . (20) 

\^Z,,...ZiJ i=i '=?--? 

This formula shows the general case with n-out going particles. The function Gi^,,,i^^ is the 
full n-point Green function where we now express it in terms of the amputated Green 
functions and only two meson propagators to suite for the current article as follows 

Summation over index i runs over two possibilities of being a neutral pion or eta. 4 and 3 
are the indices referring to the Kaon and the neutral pion respectively. Q^i is the amputated 
Green function that contains both on-shell and off-shell Feynman diagrams. The two-point 



functions are expanded near the physical poles as 

Gn{p' ^ m^ ) = -^-^^ . (22) 

The function Zi is called the wavefunction renormalization factor. The expansion of the 
off-diagonal two-point functions around the physical poles is somewhat more involved but 
can be done in terms of the one-particle irreducible two-point functions 'n.ij{iv?) and the 
mass differences needed in the propagator of the i-leg in Fig. [T] as explained in [TB] . We 
now expand all quantities to the required chiral order and use the fact that we have exactly 
diagonalized the lowest Lagrangian to obtain for the full amplitudes to order p^: 

•A.43,extv = '^43,extv + ^43,extv + •^A3,extv + " " ' ; (^<j) 

'^43,exto — yi3,extv ! (,^4) 

1 ^(4) , 1^{4)\ ^(2) ngg (3) 



,(4) 
'A3,extv 



.(4) _ ^(4) /I (4) 1 {4)\^(2) ^h8{'i) M2) (^^^ 

/I (6) _ Md) 1 fr^{6) , 7{6)\ M2) 1 fy{4) y(4)\ M 

-^43,6x41; — yi3,extv r, \"^33 "r ^44 ^ ^43,60;* r, \^^33 "r ^44 ^ i^4 

' ,4, n38(3)w ,„ 



,{2) 
'43,exti) 



n38(3)(^) ^(4) , Thsi5Y^M2) , n38(3)Wn88(3)W ^(2) 



n38(3)(^)\ ^(2) 1 /^^(4)n 38(3)W (4)n38(3)(^\ (2) 

Ami '' Ami i ^^^''^^^^ 




(26) 

The Z and 11 factors have been valuated earlier fWi [T7] and we thus need to evaluate the 
various Q amputated amplitudes. 

5 Analytical results and ratios of form-factors 

To do the calculation, we need to calculate the tree level diagrams of Figj2l the one- and 
two-loop diagrams of Fig. |3] and the two-loop diagrams with overlapping divergences of 
Fig. m with isospin breaking kept in the masses and vertices. These amplitudes should 
then be put together with the wave-function renormalization and mixing effects given in 
(123|) . The lowest order expressions are quite simple. The form- factors f!^'"^^ all vanish 
and the others are 

f^*^\t) = cos e + Vs sine, 

/f^^t) = 1, 
ff-^t) = 1, 

ff^\t) = cos e-V3 sine. (27) 
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Figure 2: The tree level Feynman diagrams for the Kaon transition form-factors. The 
wiggly line indicates the insertion of the vector current, a dot an order p^ vertex, a cross 
an order p^ vertex and a crossed circle an order p^ vertex. 

The NLO expressions agree with the isospin breaking ones calculated in [U [2] for the /+ 
form-factors. The isospin breaking in the /_ and /o form-factors is new. The NLO results 
are given in App. |X1 The full NNLO results are very lengthy but we have performed 
two independent calculations that are in agreement. All eight form-factors are also finite 
using the general subtractions calculated in [18]. The nonlocal divergences and the other 



quantities that can be removed using MS subtraction also cancel as they should. These 
consistency checks are described in detail in [19]. The loop integrals are computed using 
the methods described in [201 1^ . 

The main existing previous work is for Kf^^ decays. Isospin breaking to order p^ for /_|_ 
was done in [3] and the electromagnetic parts worked out in [7] to order e^p^. 

While this work was in progress, an analysis of the isospin breaking in the rare decay 
form-factors f^ '^ and f^ '^ to NLO and order e^p^ appeared [2]. They also noted 
that the relation f ll2p was satisfied but do not seem to have realized it is an immediate 
consequence of isospin. 

Isospin breaking in f^""^^ has not been discussed earlier within the context of ChPT. 

In [^ another relation valid to NLO and first order in isospin breaking was found. The 
ratio of form-factors 

is independent of momenta and can be cleanly predicted in terms of pseudo-scalar meson 
masses. The equality follows from the use of ( 1T2|) . 

Our results satisfy the relations ( iTTl) and ( IT2l) . we had to use a large number of relations 
between the various integrals to check this and obtained in this way another nontrivial 
check on our results. The NLO relation found by [3] is no longer true at NNLO. There 
are t-dependent corrections at order NNLO of alO types, pure two-loop, L[-dependent and 
Cj^-dependent ones. The relation f l28|) is also not true for the scalar form-factors Jq"^^^ {t) 
nor for f^^'^^t) already at NLO. 

We define here also two more ratios for later use, first the ratio of the two weak decay 



^That it was not valid for the Cf contributions at order p^ was also noticed in [2]. 
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Figure 3: The one- and two- loop Feynman diagrams for the Kaon transition form-factors 
without overlapping divergences. Notation as in Fig. O 
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Figure 4: The two-loop Feynman diagrams for the Kaon transition form-factors with over- 
lapping divergences. Notation as in Fig. [2l 



form-factors 

and second the ratio of the rare to the weak decay with charged pions in the final state 



^o-(i) = %^. (29) 






rK{t) = %^. (30) 



We define similarly definitions of r''(t), TQ_(t) and r^(t) for ratios of the scalar form-factors 

/o^'-^(t). 

6 Resonance estimate of the contribution from the CJ 

This contribution is the most difficult to estimate. In the isospin limit, f^'^{0) only depends 
on the combination {CI2 + C^a) ('^I: ~ fn^) [5] and its estimate is the main uncertainty 
in the chiral prediction for f^'^{0). A review can be found in [1]. The under lying reason 
for the factor [m'j^ — m^) is the AdemoUo-Gatto theorem[22]. The reasoning used there 
remains valid also in the case with isospin breaking for the form-factors that do not involve 
7r°-?7 mixing. The isospin conserving case is proportional to {rris — rh) , but including 
isospin breaking, the form-factor for the neutral weak decay is proportional to {nis — niu) 
and for the charged rare decay it is proportional to {mg — m^) . The full order p^ tree level 
contribution in these cases is again proportional to CI2 + C34, — {LD just as was found for 
the isospin conserved case in [HI [23] . 

The general method we use to estimate the C[ is of saturation by a finite number of 
resonances introduced by [SI |25] . We use the vector Lagrangian in the Proca formulation 
with parameters as determined in [T9| |5]. The scalar effect was studied in detail in [23] 
and more generally in [26]. Some problems with this procedure are discussed in [27] . 

The vector exchange contribution does not contribute to the values at t = for f^''^^ (t). 
It does however contribute strongly away from zero. The estimate we use here for the 
CJ from vector exchange is described in [5]. In particular, the same estimate is in good 



agreement with the estimate of the curvature in the pion electromagnetic form-factor which 
leads to an experimental determination of[2T] 

- 4 (Cg^g - C^q) = (0.22 ± 0.02) 10-3 (31) 

compared with a prediction of 0.26 10"'^. This is the part that estimates the contribution 
from the C^ in Fig. |6l The way we have implemented it here is via the effect on the Cf 
directly as given in |28j . 

Second, we take into account the contribution from the singlet pseudo-scalar degree of 
freedom Pi. We use the simple Lagrangian 

C^' = Id^Pid^Pi - ^M^i'Pi + idn.Pi{X-) ■ (32) 

Integrating out Pi leads to the order p'^ term with L7 and the order p^ Lagrangian 

A' = -^d^{x^)d'^{x~) with dm = 20 MeV. (33) 

The latter was rewritten in general in terms of the basis of operators of [13] in [17] . The 
result ici 

d^{X-)d^{X-) = O18 + ^Oi9 - io2o + ^021 + 2O27 + ^031 - ^032 + ^033 

-2035 + 037-^094. (34) 

The result is that the singlet Pi contributes via the order p^ constants C^ also to the isospin 
breaking in the values for f^'^^{0) but it does so only via tt^-i] mixing. The numerical 
result is 



eROTvO 



0.00065 



/f " (0) = -0.00065. (35) 
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7 Numerical results 

7.1 Input parameters 

For the masses we use the particle data book masses except for the eta where we use for 
consistency the value 547.3 MeV. The input values for the order p'^ constants LI we use 
are fit 10 of [16]. This fit used the i^e4 data from E865, and input values rris/rh = 24 and 
Pr/Ptt = 1.22. For the masses it used the physical masses. Electromagnetic corrections 



■*This was derived by the authors of [5S] but not included in the final manuscript. It also agrees with 
the expression shown by Kaiser [50]. 



to the Kaon mass were included with the estimate of the violation of Dashen's theorem of 
[31] included. An extensive discussion of this fit can be found in [29] using the older K^a 
data and working fully in the isospin limit. 

We will always quote results for the isospin conserving formulas of |i21j where the kaon 
mass is taken to be the mass of the kaon involved and similarly, we use for the pion mass 
the mass of the particle involved in the matrix element. For the results with the formulas 
including isospin breaking, we have used for the Kaons their physical masses but for both 
charged and neutral pion the same mass, since to first order in m^ — irtd these have the 
same mass. We have always taken the mass of the final state pion involved in the matrix 
element. The reason for this choice is to always have the kinematics right in the matrix 
elements. The effect of changing the pion mass can be judged by looking at the results for 
the isospin symmetric formulae which we quote for different input Kaon and pion masses 
in Tab. [H. The order p^ constants C[ have been put to zero at the scale /x = 770 MeV 
unless otherwise noted in Sect. 17.21 In Sects. 17.3^ 17.41 and [7.51 we have put the C[ at the 
value estimated by vector and singlet pseudo-scalar exchange at /i = 700 MeV. 

The main fit 10 with Dashen's violation gave mu/rrid = 0.45 while removing the viola- 
tion of Dashen's theorem gave 771^/771^ = 0.52. The standard values without order p^ and 
without violation of Dashen's theorem gave mu/md = 0.585 [16]. These values, together 
with the input value for 772s/rh correspond to sine = 0.0143, 0.0119 and 0.00986 respec- 
tively. This can be compared with the value of 0.0106 ± 0.0008 used in [2] which used the 
input neglecting order p^ effects. Note however that the recent evaluation from 7] -^ Sir 
[T7] leads to somewhat different values. 

7.2 /f^'(O) 

Here we give the results for the form-factor values at zero. In Tab. [1] we first show the 
results for the isospin conserving formula of [5]. Here the only isospin breaking effect is the 
different kaon and pion mass used as described in Sect. 17.11 The results for the charged and 
neutral weak decay are in agreement with [5]. We have in fact checked that the formulas 
including isospin breaking numerically agree with the isospin conserving formula if the 
masses are set to the same isospin conserving masses and sine = 0. As is clear from the 
numbers in Tab. [H the isospin breaking effects from varying the masses in the loops is 
quite small. 

In contrast, we have shown the equivalent set of values for our amplitudes including 
isospin violation. It can be seen that effect is much larger for the amplitudes with a neutral 
pion in the final state. That is, as can already be seen at lowest order, pion-eta mixing is 
important for this decay. The values in Tab. |2]are with mu/7nd = 0.45 or sine = 0.01429. 

To show the variation with the input for 7nu/7nd, we show in Tab. [3] using the same 
inputs as for Tab. |2] but with mu/rnd = 0.585 or sine = 0.009857. This corresponds to 
the fit for mu/rnd without violations of Dashen's theorem and using order p^ expressions. 
Our results, except for the lowest order in f l271) . are explicitly linear in sine. The numbers 
are slightly different from the preliminary results quoted in [8]. This is due to a slightly 
different way of treating the pion masses. 
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Table 1: The different contributions to f^''^\0) using the isospin conserving amphtudes 
of [5] . We have also shown the break-up of the order p^ expressions in the pure two-loop 
part and the L[-dependent part. The part depending on the C[ is not included. 
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Table 2: The different contributions to f^'^\0) using the amplitudes including isospin 
breaking. We have also shown the break-up of the order p^ expressions in the pure two- 
loop part and the L[-dependent part. The part depending on the CJ is not included. 
We used here mu/md = 0.45 corresponding to the two-loop fit of [16] including Dashen's 
theorem violations. 





fr-' 


ff-' 


fK^.^ 


rK'^T,^ 


order p^ 


1.01702 


1.00000 


1.00000 


0.98288 


order p'^ 


-0.01931 


-0.02282 


-0.02202 


-0.02675 


order p^ 


0.00986 


0.01467 


0.01395 


0.01919 


p^ 2-loop 


0.00435 


0.01142 


0.01084 


0.01815 


p^ L[-dependent 


0.00551 


0.00325 


0.00311 


0.00104 


sum of p^, p"^ and p^ 


1.00757 


0.99186 


0.99193 


0.97532 



Table 3: The different contributions to /^*'^^(0) using the amplitudes including isospin 
breaking. We have also shown the break-up of the order p^ expressions in the pure two- 
loop part and the L[-dependent part. The part depending on the C[ is not included. We 
used here irtu/rnd = 0.585 corresponding to the one-loop fit of [T6] without violations of 
Dashen's theorem. 



Using the results of Tab. [2] we can also quote numerical results for the various ratios 
defined earlier at the point t = 0. First the ratio of charged to neutral weak decay. This is 

ro_(0) = 1.02465 + 0.00587 - 0.00711 = 1.02341 , (36) 

where we see that the order jf' contributions lower the result and essentially cancel the 
enhancement from the ratio at order p^. If we add the contribution from singlet P\ exchange 
we obtain ro- = 1.024068. However, compared to the old order p^ value, we get again an 
enhancement due to the larger value of sine obtained from the order jP fit. We showed 
the contributions to the ratio at order p^, p^ and p^. This should be compared to the 
experimental ratio as determined from the global FLAVIAnet fit |32] 

ro-exp = 1 + i^su{2) = 1.0284 ± 0.0040 . (37) 

As we see, we obtain a reasonable agreement. 

We can also look at the double ratio r from (fT2|) . Our formulas satisfy it exactly. The 
main numerical source of the difference at higher orders results from the fact that we used 
a different pion mass in the denominator and the numerator. The result is 

r = 0.99918 - 0.00161 + 0.00085 = 0.99842 , (38) 

where a fairly sizable cancellation happens between the order p'^ and order p^ contributions. 
We again showed the contributions to the ratio at order p^, p^ and p^. 

The final ratio, of weak to rare decays with a charged pion in the final state is 

tk = 1.00000 + 0.00097 - 0.00080 = 1.00017 . (39) 

The three numbers in the middle part are once more the contributions to the ratio at order 
p^, p'^ and p^. Once more, we see a significant cancellation between the order p^ and p^ 
contributions. 

7.3 /f^'(i) 

In this subsection we show the results as a function of t for the the f^''^^ form-factors. We 
first show the case for the neutral weak decay in Figs. |5] and [61 Fig. |5] shows the result to 
lowest order, NLO and NNLO. It can be seen that there is a nice convergence in the entire 
region shown. 

We show the various subparts of the order p^ contribution in Fig. |6l The contributions 
shown are the two-loop contribution, the part dependent on the order p^ LECs L[ as well 
as the part that depends on the order p^ LECs C[. 

The results shown so far for f^ "^ are essentially the same as those in the isospin limit 
of [5] . We have included isospin breaking but it is a rather small effect for this form- factor. 
Rather than showing similar plots for the other three form-factors we show here the ratios 
as a function of t. First we show the variation of the full ratio r as a function of t. The ratio 
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Figure 5: The form- factor f^ ^ (t) as a function of t. Shown are the lowest order (p^ 
NLO (p^) and NNLO result (p^). Isospin breaking is included. 



r is somewhat more different from one than naively expected since we included different 
pion masses. The ratio ro_ defined in ( 12^ is shown as a function of t in Fig. |H1 It was 
shown in |1] that at NLO this ratio is independent of t. We have checked that this is no 
longer true at NNLO as is clearly visible in the figure. However, there is clearly no sign of 
an anomalously large isospin breaking effect in this ratio. We also show the similar ratio 
for the charged rare to neutral weak decay, r^ as defined in ( 130|) in Fig. [91 



7.4 ff^'(t) 

In this subsection we show the results as a function of t for the the f^'^^ form-factors. We 
first show the case for the neutral weak decay in Figs. [TU] and [TTJ Fig. [TD] shows the result 
to lowest order, NLO and NNLO. It can be seen that there is a nice convergence in the 
entire region shown. 

We show the various subparts of the order p^ contribution in Fig. [11] The contributions 
shown are the two-loop contribution, the part dependent on the order p^ LECs L[ as well 
as the part that depends on the order p^ LECs CJ. The latter is essentially zero here since 
the vector exchange contribution to the scalar form-factor vanishes to the order considered 
here and the singlet pseudo-scalar doesn't contribute either. A scalar exchange would 
contribute but we have not included such an estimate here. The curvature visible is here 
mainly coming from the loops. 

The results shown so far for f^ '^ are essentially the same as those in the isospin limit 
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Figure 6: The form-factor f^ ^ (t) as a function of t. Shown are the full order p^ contri- 
bution and its three constituent parts, the pure two-loop contribution, the L[-dependent 
part and the C[-dependent part. The contribution to the quadratic slope comes mainly 
from the C^ dependent part but that is fixed from the pion electromagnetic form-factor 
pT] . Isospin breaking is included. 
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Figure 7: The ratio r as defined in ( !T2|) as a function of t. Both the deviation from 1 and 
the t dependence are effects of higher order in isospin breaking. 
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Figure 8: The ratio ro_ as defined in f l29p as a function of t. This is the ratio of the 
charged to neutral weak decay. The t dependence for the NLO result is higher order in 
isospin breaking but is first order at NNLO. 
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Figure 9: The ratio tk as defined in ( 1301) as a function of t. This is the ratio of the charged 
rare to neutral weak decay. 
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Figure 10: The form-factor Jq '^ (t) as a function of t. Shown are the lowest order [j? 
NLO (p"^) and NNLO result (p^). Isospin breaking is included. 
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Figure 11: The form- factor j^ '" (t) as a function of t. Shown are the full order p^ contri- 
bution and its three constituent parts, the pure two-loop contribution, the /.[-dependent 
part and the C[-dependent part. Isospin breaking is included. 



of [5] . We have included isospin breaking but it is a rather small effect for this form- factor. 
Rather than showing similar plots for the other three form-factors we show here the ratios 
as a function of t. First we show the variation of the full ratio r° as a function of t. The 
ratio r^ is somewhat larger than naively expected since we included different pion masses 
The ratio rQ_ defined in ( 129|) but for the scalar form-factor is shown as a function of t in 
Fig. [131 This ratio can be t-dependent already at NLO which is clearly visible. However, 
there is no sign of an anomalously large isospin breaking effect in this ratio. We also 
show the similar ratio for the charged rare to neutral weak decay, r^ as defined in fl30p in 
Fig. [3 The scalar form-factors are not needed for the weak decays to an electron or the 
rare decays to a neutrino-antineutrino pair. They do contribute to the weak decays to a 
muon and the rare decays with a muon-anti-muon pair via a the axial current couplings to 
the latter from the short- distance contributions. 



7.5 Callan-Treiman point 

The Callan-Treiman relation [33] states that the scalar form-factor at t 



m 



^--m? satisfies 



K 



/o( 



m^ — m^ 



+ 0{mu,md). 



(40) 



This relation is derived using current algebra in the up-down-sector and should thus have 
rather small corrections of order m^. The relation is exact when the up and down quark 
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Figure 12: The ratio r^ as defined in flT2|) bit for tlie scalar form-factor as a function of 
t. Both the deviation from 1 and the t dependence are effects of higher order in isospin 
breaking. 



masses are zero. The correction at NLO was worked out in [1] and found to be 

Act ^ h (< -ml)-^ = -3.5 10-^0] . (41) 



The correction in the isospin hmit at NNLO was never presented in [5j. We have calculated 
this and also for the four different amplitudes at order p'^ but for the C[ estimates we have 
used the vector and pseudo-scalar singlet contributions as described in Sect. El 

The inputs we used produce Fk/F^, = 1.22, which is what we have subtracted from the 
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in the numbers quoted below. The isospin symmetric expression with 



and mi 



mz 



gives 



A 



CT 



-6.210" 



(42) 



As we see there is a substantial NNLO correction. Note that we did not include the 
contributions from nonzero C[2, Cu i^ ^^^^ expression. These read 



16 



^CT\cr = TrA^Ci, + C;,)ml{ml 



m. 



(43) 



It is clear that this satisfies the Callan-Treiman theorem. Notice that it is the same 
combination of order p^ LECs that shows up in the scalar slope when the part via F^/F^^ 
is subtracted as in Eq. (5.2) of [5]. 
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Figure 13: The ratio rQ_ as defined in (129|) but for tfie scalar form-factor as a function of 
t. This is the ratio of the charged to neutral weak decay. The t dependence is first order 
in isospin breaking both at NLO and NNLO. 
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Figure 14: The ratio r^ as defined in ( 130|) but for the scalar form-factor as a function of t. 
This is the ratio of the charged rare to neutral weak decay. 



For the expression including isospin breaking we simply present the numerical results 
directly 

A^;-° = 15.110-3, 

A^°-" = -5.610-^ 

A^+-+ = -9.410-^ 

A^""" = -26.4 10-^ (44) 

Recently, Leutwyler discussed the experimental measurements and the extrapolation to 
the Callan-Treiman point [31] . The results we obtained are clearly not sufficient to explain 
the large value of A^r = -0.071 ± 0.014^^48 ± 0.002^^^0 ± O.OOSe^t observed by NA48 in 
the charged weak decay [35] but are in reasonable agreement withe one observed by KLOE 
for the neutral weak decay [36] . 

8 Conclusions 

In this paper we have calculated all the vector form-factors of Kaon to Pion transitions to 
ffist order in the quark mass difference m„ — rrid and to NNLO in ChPT. We have thus 
calculated the eight different form- factors defined in Eqs. ([2ll5]) to order p^{mu — md). This 
complements the earlier calculations to order p'^{mu — w^d) done for the /+ form- factors in 
[1] and [2] and to order p^ in the isospin limit for the vector and scalar form-factor. 

What we find in all cases is the the NNLO results diminish the effects of isospin breaking 
but due to the change in mu/md from a NLO to a NNLO fit the total effect is to increase 
the isospin breaking in the form-factors. This goes some way towards reconciling the 
determinations of Vus from the charged and neutral weak K^y, decays but does not explain 
the full difference. We have also calculated isospin breaking in all the scalar form-factors. 
Here again, the effects are sizable but not unexpectedly large. In particular, they are 
not large enough to explain the discrepancy with the Callan-Treiman point observed by 
NA48[35]. 

We have presented numerical results for the values at t = and for the t-dependence 
as well as for various ratios of the form-factors. In particular, we have shown that the 
relations ( ITTl) and ( IT2l) are valid to all orders in ChPT and to ffist order in ttIu — rrid. We 
presented numerical results for the ratios r, ro_ and tk as well as for their equivalents for 
the scalar form-factors. 
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A The order p^ expression 

In this appendix we explicitly write out our order p^ results. We have checked that they 
agree with the published results for the /+ form- factors of [1] and [2]. They also satsify 
the relation ( [T2|) when the integrals are expanded to obtain a common Kaon mass, but we 
have quoted all eight formf actors here since by rewriting one can move things between the 
order p^ and p^. The expressions quoted here are the ones we used to define the order p^ 
part. The integrals used below are the standard one- loop integrals defined in many places, 
see e.g. [20] . 
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+ Si(m2o, m^o, t) (+1/12 1 + 1/12 m^o - 13/12 m^o) 
+ Si(m^o,mJ,t)(-l/4t + 23/12m|^o- 11/12 m^o) 
+B2i{ml+,mj^+,t){-t- m|-o + m^o) 
+ S2i(m^o,m|^o,t) (+l/6t- l/2m^o + 1/2 m^o) 
+ 52i(m^o,m5,t)(-l/2t-3/2m^o + 3/2m2o) 

-522(m^+,m^+,t) + l/6B22imlo,m%o,t) - l/2 522(m^o,mJ, t)) . (45) 



.6 



B The order p^ LECs dependent part 

In this appendix we write out explicitly the part dependent on the order p^ LECs C[. We 
use here a notation which uses the property (ITU|1 . 






ff^'it) 



C" 



c 



c 



c: 



We also use the notation 






Sine 



V3 ( m^o - m 



ff{t) 



sine 



^(//(t) ^//^(t) ^(^2^^_ 



2 2,2 2 



T?2" 



/f(t) 



(46) 



(47) 



The pion mass we have used generically since they are the same to the order of our calcu- 
lation. 

The CI dependence is now given by 



m 



-e (-4 CI, + 4 C;,) + ml t (-4 C[2 - 16 CI, - 4 Q3 - 4 Q^ - 2 q; 



90/1 



+m: t (-12 CI2 - 32 C[3 - 4 Q3 - 8 Q4 - 4 Q^ - 6 C^.^) + < (-2 ^[2 - 2 ^3^,) 
+m^ m^ (4 C{2 + 4 Cg^J + m^ (-2 ^[3 - 2 Cg^J , 



fB^t) = +t2(_i2Q3 + i2C9g+m^(-4C[2-48C[3-4Q3-12Q4-2Qo) 
+ml t (-44 CI, - 96 ^[3 - 20 C^, - 24 C^, - 12 C^, - 22 C;,) 
+ml (2 CI, + 16 C[4 + 16 CI7 + 48 ^[3 - 14 C^^ - 24 C^g) 
+ml ml (-4 C[2 - 32 ^[4 - 32 ^[7 - 96 Cl^ + 28 ^3^4 + 48 ^3^5) 
+mt (2 CI, + 16 C[4 + 16 Clr + 48 C^^ - 14 ^3^4 - 24 ^3^5) , 

/+W = +m^ (96 C[9 + 64^2^0 + 64^3^1 + 64^3^2 + 128 ^3^ 

+ml ml (-32 ^[4 - 32 ^[7 - 96 Cl^) + m^ m^ (64 Cl^ + 64 C[7 + 192 C{^ 

-288 C[9 - 192 CIq - 192 ^3"! - 192 CI, - 384 C3"3) 

+ml (-32 C[4 - 32 C[j - 96 C[g + 192 ^[g + 128 C^q + 128 CJ^ + 128 C;, 

+25QC;,), 

f^{t) = +mlt{8Cl,-8C^,, + 8C^,, + 4C^,o)+mlt{-8Cl, + 8C^,,-8C;,-4C 



r \ 
90 J 



+ml (8 CI, + 8 CI,) + ml ml (-16 ^[^ - 16 ^3^4) + m^ (8 CI, + 8 ^3^4) , 

f^{t) = +mlt{-4Cl, + 2C;,-2C;,)+mlt{ACl,-2C;, + 2C^,o) 

+mt{6Cl, + 8Cl, + ACl, + ACl, + 2C;, + 2C;, + 2C;, + C;,) 
+mlmlil2C{, + 8C{, + AC{, + 8C{, + AC!,, + 2C^,, + 2C^,, + 2C^,,) 
+mt (-18 CI, - 16 CI, - 4 C[4 - 8 CI, - 8 CI, - 6 C^4 - 2 Q3 - 4 Q4 

~2 Cg5 — 3 Cgg) , 

/f(t) = +m^ (-4 C[2 + 2 Qs - 2 eg + m^t (4^2-2^8 + 2^0) 

+2 Cg4 — C90) 

+ml ml (36 C[2 + 8 ^[3 + 16 ^[4 + 4 ^[5 + 72 C^j + 96 C[8 + 44 C3"4 + 48 C;, 

+8 Cgs + 2 ^64 + 2 Cgs + 6 C90) 

+mt (-30 CI, - 16 C[3 - 12 C^, - 8 ^[5 - 40 ^[7 - 48 ^[g - 26 ^3^4 - 24 C;, 

f^{t) = 0, 

f^{t) = +mlt{8Cl,-4C;, + AC;,)+mlt{-8Cl, + 4C;,-4C;,) 

+mt (-24 CI, - 16 C[3 - 8 CI, - 16 ^[7 - 8 ^3^4 - 4 C^, - 4 C^^^ - 4 Cg^^) 

+m: m^ (-16 C{, - 16 ^[4 - 8 CI, + 16 ^[7 - 8 Q3 - 4 Q4 + 4 Q^) 

+mt (24 CI, + 32 ^[3 + 16 CI, + 16 CI, + 8 ^3^4 + 8 Q3 + 8 Q4 + 4 Cg^o) . (48) 

References 

[1] V. Cirigliano, PoS KAON (2007) 007. 

[2] F. Mescia and C. Smith. EFXw: 0705 .20251 [hep-ph]. 

[3] J. A. Cronin, Phys. Rev. 161 (1967) 1483. 



[4] 
[5] 
[6] 
[7] 

[8 
[9 

[lo: 
[11 

[12 
[13 
[14 

[15 
[16 

[17 
[18 

[19 

[20 

[21 
[22; 



J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 517. 

J. Bijnens and P. Talavera, Nucl. Phys. B 669 (2003) 341 |arXiv:hep-ph/0303103| . 



P. Post and K. Schilcher, Eur. Phys. J. C 25 (2002) 427 tarXiv:hep-ph/0112352| . 

V. Cirighano, M. Knecht, H. Neufeld, H. Rupertsberger and P. Talavera, Eur. Phys. 
J. C 23 (2002) 121 |arXiv:hep-ph/01 101531 . 

J. Bijnens, PoS KAON (2007) 027 [ arXiv:0707.04T9l [hep-ph]]. 

V. Bernard, M. Oertel, E. Passemar and J. Stern. 'arX iv:0707.4194l [hep-ph]: 
Phys. Lett. B 638 (2006) 480 | arXiv:hep-ph/0603202] . 

M. Jamin, J. A. Oiler and A. Pich, JHEP 0402 (2004) 047 [arXiv:hep-ph/ 0401080l . 

S. Weinberg, Physica A 96 (1979) 327. 

J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

J. Bijnens, G. Colangelo and G. Ecker, JHEP 9902 (1999) 020 



arXiv:hep-ph/9902437i . 
J. Bijnens, Prog. Part. Nucl. Phys. 58 (2007) 521 [arXiv:hep-ph/0604043| . 



G. A moros, J. Bi jnens and P. Talavera, Nucl. Phys. B 602 (2001) 87 
(aBGv:hep-ph/0101127j . 

J. Bijnens and K. Ghorbani. larXiv:0709.02"30l [hep-ph]. 

J. Bijnens, G. Cola ngelo and G. Ecker, Annals Phys. 280 (2000) 100 
arXiv:hep-piV9907333| . 



J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and M. E. Sainio, Nucl. Phys. B 508 
(1997) 263 [Erratum-ibid. B 517 (1998) 639] [a rXiv:hep-ph/9707291] . 

G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 568 (2000) 319 
[arXiv:hep-ph/9907264] . 

J. Bijnens and P. Talavera, JHEP 0203, 046 (2002) |arXiv:hep-ph/0203049| . 

M. AdemoUo and R. Gatto, Phys. Rev. Lett. 13 (1964) 264; 
R. E. Behrends and A. Sirlin, Phys. Rev. Lett. 4, 186 (I960). 



[23] V. Cirigliano, G. Ecker, M. EidemuUer, R. Kaiser, A. Pich and J. Portoles, JHEP 
0504, 006 (2005) [arXiv:hep-ph/0503108| . 



[24] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[25] G. Ecker, J. Gasser, H. Leutwyler, A. Pich and E. de Rafael, Phys. Lett. B 223 (1989) 
425. 

[26] V. Cirigliano, G. Ecker, M. Eideniuller, R. Kaiser, A. Pich and J. Portoles, Nucl. Phys. 
B 753 (2006) 139 [a rXiv:hep-ph/0603205] . 

[27] J. Bijnens, E. Gamiz, E. Lipartia and J. Prades, JHEP 0304 (2003) 055 
(arXiv:hep-ph/0304222j . 

[28] K. Kampf, J. Novo tny and J. Trnka, Eur. Phys. J. C 50 (2007) 385 
|arXiv:hep-ph/060805li . 

[29] G. Amoros, J. Bijn ens and P. Talavera, Nu cl. Phys. B 585 (2000) 293 [Erratum-ibid. 
B 598 (2001) 665] |arXiv:hep-ph /0003258|. 

[30] R. Kaiser, talk presented at EuroFlavour06, 2-4 November 2006, Barcelona, Spain. 



[31] J. Bijnens and J. Prades, Nucl. Phys. B 490 (1997) 239 |a rXiv:hep-ph/9610360| . 

[32] M. Palutan, PoS KAON (2007) 020 

[33] G. G. Gallan and S. B. Treiman, Phys. Rev. Lett. 16 (1966) 153. 

[34] H. Leutwvler. larXiv:0706":3T38l [hep-ph]. 



[35] A. Lai et al. [NA48 Collaboration] , |arXiv:hep-ex/0703002 



[36] F. Ambrosino et al. [KLOE Collaboration]. iarXiv:0707.4631 [hep-ex]. 



